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1 Introduction 



Although the parton model has been very successful in providing a qualitative 
description of e + e~ annihilation to hadrons, deep inelastic lepton-hadron scat- 
tering (DIS), and Drell-Yan lepton pair production, a more detailed description 
requires higher-order QCD corrections. These have been known for many years 
to first order in a s , and more recently to second order, but are generally presented 
in an inclusive form. However, in many instances, a more exclusive treatment is 
needed, for example to describe features of the hadronic final state or correlations 
between the hadronic and leptonic parts of the process. This is particularly true 
of Monte Carlo event generators, where the aim is to provide a fully exclusive 
description of the process, event by event. 

In this paper we present the first order tree-level corrections to quark scatter- 
ing and annihilation processes in a form in which the QCD and electroweak parts 
exactly factorise. This makes exclusive event features and correlations particu- 
larly transparent. We also discuss how these could be used to provide a simple 
Monte Carlo treatment of the processes. We essentially follow the method of 
Ref. |1| for e + e~ annihilation, taking over most of the same notation. 

In anticipation of the main applications of our results, we use the language 
of the DIS and Drell-Yan processes, but in fact our method is applicable to any 
electroweak process in which the lowest-order diagrams contain a single quark 
line attached to a single gauge boson. It should also be stressed that the method 
treats the quarks as massless throughout. 

The paper is set out as follows: In the remainder of this section, we recap the 
important ingredients of [|l[] . In section ^| we follow the same method to derive the 
equivalent expression for the QCD Compton part of the first-order correction to 
DIS (QCDC, q£ — > qg£). In section |3| we do the same for the boson gluon fusion 
part (BGF, g£ — > qq£). In section |] as an example of the simplicity of our form, 
we discuss lepton-hadron correlations in DIS, which have been proposed as an 
important test of QCD[[J,[J. In section ^ we derive and discuss the equivalent 
results for the Drell-Yan process. Finally in section [5] we give a summary. 

The tree- level Feynman diagrams for e + e~ annihilation to hadrons at 0(a s ) 
are shown in Fig. [I]. It is well known that only the sum of the two diagrams is 
gauge invariant, and that they can never be separated. However, the first impor- 
tant step of |J is to use the explicit gauge choice introduced by the CALKUL 
collaboration ^J], in which the diagrams only contribute to physically distinct pro- 
cesses, with all the interference absorbed into the form of the polarisation tensor. 
In this gauge, only the first diagram of Fig. [1] contributes to final states in which 
the gluon has the same helicity as the quark, and the second diagram when it has 
the same helicity as the antiquark. The polarisation tensor contains collinear- 
divergent pieces, so that both diagrams still give collinear divergences in both 
directions. 

The next important step is to introduce a pair of massless four-vectors r\ % i 
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Figure 1: Tree-level Feynman diagrams for e + e~ — > hadrons at 0(a s ), in which 
the lepton side is represented by an arbitrary current, J M , and the boson-quark 
coupling by uj^ . 

that are parallel to px i2 but with energy Q/2, the energy of the partons in the 
parton model. These allow each matrix element to be written as the product 
of two parts, one corresponding to parton-model production, the other to the 
emission of the gluon without any knowledge of J M or uA The final result is then 

das = ^dT 2 dXl ^ 2 f {xl\M 2 (r u q - r,)| 2 + x\\M 2 {q - r 2 , r 2 )| 2 }, 
27r (1 — Xi)(l — x 2 ) 2n 1 > 

(1) 

where Xi = 2pi-q / q-q are the energy fractions in the cmf of q, dT 2 is an element 
of the parton-model phase-space, and <fi is the azimuthal angle of the gluon. Both 
dT 2 and 4> have different interpretations for the two parts of the matrix element: 
in the first they refer to the phase-space to produce a quark with momentum 
r\ in the parton model, and the azimuthal angle around r±, while in the second 
they refer to an antiquark with momentum r 2 . With the exception of neglecting 
quark masses, this expression is exact for arbitrary currents of vector or axial- 
vector type, so applies equally well to the W decays in WW pair production for 
example. It is also valid at the helicity level, so if polarisation effects are retained 
in the lowest-order matrix element, they are correctly treated at 0(a s ). 

As a Monte Carlo prescription, (|IJ) has a simple interpretation: first parton- 
model events are generated according to the exact matrix element, then gluon 
emission is generated according to 

1 da C F a s x\ + x\ 



a dxidx 2 2n (1 — Xi)(l — x 2 ) ' 



(2) 



and finally parton 1 or 2 is chosen with relative probability x\ and x\ to retain 
its parton-model direction, with the hadron plane being rotated uniformly in 
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azimuth about it. The quarks retain their lowest-order polarisation states, with 
the gluon inheriting the same state as the parton that was chosen. 

This procedure produces the exact distribution of three-parton final states, 
with no knowledge of the mechanism that produced them, provided an exact 
generator of two-parton final states is available. 

2 QCD Compton Process in DIS 
2.1 The Matrix Element 

The three-parton matrix element for a quark of helicity A and a gluon of helicity 
p is given by 

where the momentum of the incoming(outgoing) quark is P\{p 2 ), and the gluon 
is k. All the partons are considered massless. We use the explicit formulation of 
the gluon polarisation tensor e p (k) introduced by the CALKUL collaboration, 



^{k) = iV(i(l + P7 5 W 2 ^i-^i^(l + P7 5 )), (4) 

_i 

N = [4pi-k p 2 -k p\-p 2 ] 2 • (5) 

After a little manipulation, we obtain 

Mt = M 3 (X, A) = -ig a Nu x (p 2 )u ,i W 1 -fl& l u x (p 1 )J li (q), (6) 

M 3 =M 3 (X,-X) = -igsNu x (p2)AW 2 + V)w' i u x (p 1 )J ll (q). (7) 



The special vectors we introduce are defined covariantly by the momentum frac- 
tions of the partons, 



just as in e + e annihilation. In general X\ < —1, x 2 < 1, where the equalities 
apply in the parton-model scattering case. 



n = -Pi/xi, (9) 

r 2 = n + q = p 2 + k-p 1 -p 1 /x 1 , (10) 

r 2 = p 2 /x 2 , (11) 

fi = r 2 -q=p 1 -k-p 2 +p 2 /x 2 . (12) 



It can be seen that r±, r 2 are the momenta that the incoming and outgoing quarks 
would have in 0{\) scattering at the same y B and Q 2 . fi,r 2 have the same 



2p 2 -k 



-2p r k 



u x (pi)J^(q), (3) 
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component parallel to the current direction, but also have a component transverse 
to it, such that r 2 is parallel to p 2 . As we discuss in more detail later, these are 
the momenta that the incoming and outgoing quarks would have in the parton 
model with an 'intrinsic' pt- 

We can use these vectors to rewrite the matrix elements using the following 
results, 

u(ap) = \/au(p), (13) 
Wi-flti = fijfc (14) 
= u x (ri)u x (ri)^2, (15) 
tfiifa + ft) = tfiU\(f 2 )u\(f 2 ), (16) 

where = denotes "equal except for an overall complex phase" , and all the results 
rely on the assumed masslessness of pi )2 . Inserting these into the matrix elements 
we obtain 





= g s N^fx2~ u x (r 2 )u^u x {r 1 )u x {r 1 )p ! 2 u x {p l )J^q) 


(17) 




= C + M 2 (r 1 ,r 2 ), 


(18) 


C + 


= gsNy/x^uxif^uxipi), 


(19) 


M z 


= C~M 2 (r 1 ,f 2 ), 


(20) 


C~ 


= 9sNy/xlux(p2)¥llU\(f2), 


(21) 



where Ai 2 (qi,q 2 ) = u x {q 2 )u^u x {qi)J^{q) is the parton-model matrix element for 
the current to scatter an incoming quark q±, to an outgoing quark q 2 . Finally we 
can calculate the \C ± \ 2 explicitly to give 

|C-| 2 = ^ —x\. (23) 



2.2 Kinematics and Phase-Space 

We parametrise the parton model scattering by x B , Q 2 and $, the lab-frame 
azimuth of the lepton scattering, which is uniform for leptons with no trans- 
verse polarisation. In addition, we parametrise the hadron-plane momenta by 
X\ >2 and 0, the azimuth of the outgoing quark around the boson-hadron axis in 
their cmf (or, equivalently, the Breit frame). To be precise, when discussing the 
specific case of DIS, = when the outgoing lepton and quark are parallel. The 
phase-space limits are 

— < x x < -1, (24) 
l+x 1 < x 2 < 1. (25) 
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The cross-section at 0(1) is then 



do 2 = ^^/(*B,QW^B^£|.M 2 r (26) 
= dT 2 H\M 2 \ 2 , (27) 



and at 0(a s 



= 128^ ^/(-^ 1 ,Q 2 )<iQ 2 <ix B ^^ 2 ^Q 2 E|A1 3 | 2 (28) 

= I^^ " r,I '»" I S Q 2 Z\M 3 f (29) 

4(27r) J x B f(x B ,Q 2 ) -xi 2tt 

C f q; s — x B Xif(— x B Xi, Q 2 ) dxxdx, d(f> 

-di 2 



2tt " x B f(x B ,Q 2 ) -x\(-l - xi)(l - x 2 ) 2-k 
{xl\M 2 (n,f 2 )\ 2 +x 2 2 \M 2 (f 1 ,r 2 )\ 2 } . (30) 

It is immediately clear that this is closely related to the e + e~ formula (|ip, and is 
also an exact factorisation of the QCD and electroweak parts of the process. 

However there is one important difference from the e + e~ case: while the 
parton-model configurations described by r\ and r 2 are related by a rotation 
in e + e~, they are related by a transverse boost in the scattering process. This 
means that the simple Monte Carlo prescription is not possible because although 
it seems natural to generate lowest order e + e~ annihilation events with all possible 
orientations, it is not so natural to generate scattering events with all possible 
intrinsic transverse momenta. 

It is worth noting that x 2 is not constrained to be positive, so the vector r 2 can 
have a negative energy-component. If one forgets this fact, and simply evaluates 
the resulting matrix elements, then (^) works correctly, but to provide a physical 
interpretation, one should first use CPT- invariance to rewrite the matrix element 
to scatter a quark from f\ to r 2 as the matrix element to scatter an antiquark 
from —r 2 to — f\. 

In generating events, it is more convenient to introduce different phase-space 
variables, x p = —1/xi and z p = = 1 + (1 — x 2 )/xx, which have independent 
phase-space limits 

x B < x p < 1, (31) 
< z p < 1. (32) 



Inverting the relations we get 



xt = — -, (33) 

Xp 

x 2 = 1-— ^. (34) 

Xp 
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We also make use of x±, the transverse momentum, measured in units of Q/2, 

x 2 = — Xp ^^ — Zp ^ Zp (35) 

Xp 

to rewrite the expressions for |C ± | 2 , 



|/^»+|2 87TQ: S / 2/ 2 , 2\1 ^2 / Qfi \ 

|C 1 " (1 - x p ){\ - z p )Q* W* 2 + X±) f xJTxJ' {db} 

^ ? = {l-x p ){l-z p )Q* {l} - (37) 

The reason for writing \C + \ 2 in this form is that as x 2 tends to zero, the energy 
of r 2 tends to infinity, and the matrix element \M. 2 {fi,r 2 )\ 2 in general diverges 
as x\fx\, so the product 

^Vl-^(n,r 2 )| 2 (38) 

remains finite. The expressions in curly brackets in \C ± \ 2 both lie between zero 
and one throughout the physical phase-space. 
The cross-section is then 

C F a s i^f^iQ 2 ) dx p dz p d(f) f 2 

(^(^ + 4))-rf-rlM2(f 1 ,r 2 )| 2 l (39) 

The product dT 2 \A4 2 (ri, r 2 )\ 2 is exactly the lowest order differential cross-section, 
so corresponding events are provided by a parton-model event generator, but 
\M- 2 {fi,r 2 )\ 2 requires us to reevaluate the lowest order matrix element for the 
new momenta. 

The Monte Carlo procedure is then as follows: 

1. Generate a parton model event according to the exact 0(1) matrix element. 

2. Generate x p and z p values according to dx p dz p /((l — x p )(l — z p )). 

3. Generate a value uniformly, and construct the corresponding momenta 
Pi, P2, k and n )2 . 



4. Calculate the ratio 



R 2 s ^ ^ftjf . (40) 
x 2 2 + x\ \M 2 (r u r 2 )\ 2 
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5. Calculate a weight factor 



^^fj^m^ + X * i4 + X ^' <41) 

6. Keep the event with probability proportional to w. 

This produces the exact distributions of x p , z p and 0, for the given lowest-order 
phase-space point. The average value of the weight factor is the total 0(a s ) 
correction to the cross-section (at tree-level, according to the chosen cutoff). 

This is as far as one can take the exact calculation for a general quark scat- 
tering process, but if we now specialise to the DIS case, we can proceed further. 
The ratio R 2 can be quite generally written 

cos 2 6 2 + A cos 6 2 (l — \Jl 2 — 1 sin 9 2 cos (f^j + (l — \Jl 2 — 1 sin 9 2 cos 

R2= 1 + Al + l 2 ' 

(42) 

where cos# 2 = x 2 / ' \jx\ + x \i sin6> 2 = x±/ ' \jx\ + x 2 ± , and I = 2/y B — 1, are all 
Lorentz-invariant quantities that have simple interpretations in the Breit frame: 
62 is the angle between P2 and the exchanged boson direction, and Ql/2 is the 
energy of the incoming lepton. A is related to the couplings of the fermions to 
the exchanged bosons, 

(Cy,e + C\ ( ){C Vq + C 2 Aq ) 

For pure photon exchange A^ = 0, for charged- current interactions A cc = 2, and 
the full neutral-current case including 7-Z interference can be calculated from 

cv„ = e.+(W2-e,rf«^|^-, (44) 

Q 2 

Ca,i = h,i /2 lna o \ . a (45) 

(Q 2 + m 2 .) sin 9 W cos 6 W 

Polarisation effects can also be incorporated by using appropriate couplings in A. 
Using these expressions we can take the azimuthal average of R 2 analytically and 
generate the distribution exactly, improving the Monte Carlo weight distribu- 
tion. It also allows a direct interpretation of the azimuthal effects, as we discuss 
in a later section. 



3 Boson Gluon Fusion 

The treatment of boson gluon fusion is mostly similar to QCD Compton, so we 
do not go into as great detail. 
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The matrix element is 
M 3 (\,p) = ig s u\(p 2 ) 



Ux(pi)Jp(q), (46) 



where the momenta of the quark and antiquark are p 2 and p 1 , and all else is as 
before. This can be simply obtained from the QCD Compton case, by crossing 



Pi 
k 



~Pi, 
-k. 



The special vectors are 

r\ 

r 2 
r 2 
fi 



ri + q = p 2 ~k + p 1 -p 1 /x 1 , 

P2/X 2 , 

r 2 - q = -p\ + k -p 2 + p 2 /x 2 . 



(47) 
(48) 

(49) 
(50) 
(51) 
(52) 



Note that X\ now applies to the outgoing antiquark, rather than the incoming 
quark. This time, neither set corresponds to lowest-order scattering — both have 
transverse components. We obtain 



|C ' " (1 - Xl )(l - x 2 )Q* X » (53) 

|C ~' 2 = (1 - Xl )™- xJQ*^ (54) 

We again parametrise the kinematics by x v and z p , where z p still refers to the 
outgoing quark, ie. particle 2, and x p refers to the incoming particle, which is this 
time the gluon. We then obtain 



x 2 



1 - 



x p 
1 - z„ 



x r . 



This gives 

|C+| 2 = 

|C-| 2 = 
The cross-section is then 



in a. 



da 



dT 2 



z p (l-z p )Qt& {xl + xl) } 
z p {\-z S p )Q^ 

p/sdf'^ 2 ) dx p dz p d(f) 



x 2 



2 i 2 ' 
X 2 + x\_ 



x\ 



x\ + x\ ' 



(55) 
(56) 

(57) 
(58) 



•'T! 



2tt x B f q (x B , Q 2 ) z p (l - z p ) 2tt 
{(xl(x 2 + xl)) R, + (x 2 p (xl + xl)) R 2 ) \M 2 ( qi ,q 2 )\ 2 , (59) 
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where q\^ 2 are the parton-model momenta, and 



A \M 2 {r 1 ,f 2 ) 

x\ + x\ \M 2 (qi,q2) 
x\ \M 2 (fi,r 2 y 2 



R 1 \ JV1 2K<1,<2)\ 

1 -2 , _2 I KAJn, ^12' 



+ |A< 2 (gi,fc)| 2 ' l ° J 

which can be rewritten for the explicit case of DIS as before. We obtain the 
identical expression for R 2 , and the same for Ri, but with 6 2 and 4> replaced by 
the corresponding Q\ and 7r — <fi. 

Up to here, the BGF cross-section has been written purely in terms of the 
lowest-order quark scattering (ie. not antiquarks). There is nothing wrong with 
this in itself, as the exact distributions of BGF events will be correctly reproduced. 
However, if we want to view this as a correction to the lowest-order process, rather 
than just a cross-section calculation, it would certainly be more desirable to treat 
quarks and antiquarks equivalently. We can do this by noting that the z p = 1 
singularity is associated with configurations that become collinear to lowest-order 
quark scattering, while z p = is associated with antiquark scattering. We use 
the separation 

— tt^-t T = — + (62) 

ZpyL Zp) Zp 1 Zp 

to rewrite the cross-section 

, \a s %fA%iQ 2 )dXpdZp d(f) 

da = ±—dT 2 p v y — 

2n x B fg{x B ,Q 2 ) I - z p 2n 

{{x 2 p{x\ + xl)) R 1 + (x 2 p{x\ + xl)) R 2 ) \M 2 {qi,q 2 )\ 2 



(\, S/s(f?>Q ) dx v dz v 



2tt x B f g (x B ,Q 2 ) z p 2tt 
{(x 2 p (xl + xl)) Rt + (x 2 p(x 2 2 + xl)) R 2 ) \M 2 ( qi ,q 2 )\ 2 . (63) 

Now we define a set of exchanged variables, 

Pi = P2, (64) 
P2 = Pi, (65) 
z p = 1 - z p . (66) 

Although the special vectors are defined in the same way, they become negated, 
because of the minus sign in ([49]), so we obtain 



da ^ as dY ^ dx p dz , 



2tt x B f q (x B ,Q 2 ) I- z p 2tt 
{ (x 2 p {x\ + x\)) R t + (x 2 p {x\ + x\)) R 2 ) \M 2 (qi, q 2 ) 



10 



_i_ i^±dr- ^ dx v dz p 



2tt x B f q (x B ,Q 2 ) 1-Zp 2ir 
(x 2 v {x\ + xl)) —^--\M 2 (-r 2: -~f 1 



+ + 4)) 3^3-|^2(-F 2 , -fi)| 2 } . (67) 

Finally, we can use the CPT-invariance of the matrix element, to write 

\M 2 (-q 2 ,- qi )\ 2 = \M 2 ( qi ,q 2 )\ 2 , (68) 

where Ai 2 (qi, q 2 ) is the matrix element to scatter an antiquark from qi to q 2 . Thus 
the two halves of the cross-section are identical, but with quarks and antiquarks 
interchanged, 

da ^-dT ^ g ^ ,( ^ dxpdzp dxp 
2tt 2 x B f q (x B ,Q 2 ) 1 - z p 2tt 

{(xl{x\ + x\)) R x + (x 2 p{xl + x\)) R 2 ) \M 2 ( qi , q 2 )\ 2 
{-ot s l^fgil^iQ 2 ) dxpdzp d(f) 



2ir x B fq(x B ,Q 2 ) l-z p 2n 
{(x 2 p(xl + x 2 ± )) R, + (x 2 p(x 2 2 + x 2 ± )) R 2 ) \M 2 ( qi ,q 2 )\ 2 . (69) 

The two halves can then be separately associated with lowest-order scattering of 
quarks and antiquarks. 

The Monte Carlo algorithm is then: 

1. Generate a parton model quark or antiquark event according to the exact 
0(1) matrix element. 

2. Generate x p and z p values according to dx p dz p /{l — z p ). 

3. Generate a <fi value uniformly, and construct the corresponding momenta 
Pi, p 2 , k and n )2 . 

4. Calculate the ratios R\ and R 2 . 

5. Calculate a weight factor 

1 2b f f£B £)2\ 

- - + + 4(4 + 4) * 2 > ■ (70) 

6. Keep the event with probability proportional to w. 

This again produces the exact distributions at the given phase-space point, and 
gives the total tree-level 0(a s ) correction there. 
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4 Lepton-Hadron Correlations in DIS 



Before discussing the azimuthal correlations that arise in first order QCD, it is 
worth recalling that a correlation also arises in the parton model, if the parton's 
intrinsic p t in a hadron is included. This is because at fixed Q 2 , different azimuths 
correspond to different parton-lepton collision energies, so even though the Fermi 
motion of partons in a hadron is azimuthally uniform, their scattering cross- 
section is not. The matrix element is identical to \M.2(fi, ?"2)| 2 evaluated earlier, 

cos 2 82 + A cos 62 (l — Vl 2 — 1 sin 62 cos 6] + (l — y/l 2 — 1 sin 82 cos 6] 

l^l 2 K " TTaT+p L 

(71) 

It is this partonic matrix element that determines the azimuthal asymmetry 0, 
and not the corresponding partonic cross-section (~ |A^| 2 /s, where s is the 
parton-lepton invariant mass). 82 can be directly calculated since in the Breit 
frame, the incoming parton has longitudinal momentum Q/2 and transverse mo- 
mentum p t , so we have 

tan# 2 = 2p t /Q. (72) 

The size of the correlation is usually parametrised by the average values of cos 6 
and cos 26, 



-A cos 8 2 VP^T sin 8 2 - 21 s/F^l sin 8 2 
(COS0) ~ 2(cos 2 82 + Acos8 2 l + l 2 + \(l 2 - 1) sin2 02) ' 

|(l 2 -l)sin 2 8 2 

( 0) ~ 2(cos 2 82 + Acos8 2 l + l 2 + \{l 2 - 1) sin 2 8 2 ) ' ( } 

It can be seen that both effects are maximised by working at small Q 2 , so that 
82 is maximised (though it is still small for accessible values of Q 2 ). It is also 
increased somewhat by working at small y B (ie. / ^> 1). Taking the limit pt -C Q 
and using the definition of /, we obtain the usual expressions ||, 

(COS0) = 1 + (1-,b) 2 + ^b(2-, b ) (-Q-J' (75) 
, o ,\ 2(1 -y B ) (p 2 \ 

In the small-^B limit, these cross-sections become independent of A and, as noted 
in [|5| are even the same if the fermions are replaced by scalars. This corresponds 
to the well-known fact that the cross-section for scattering by exchange of a single 
spin-1 particle has the same high-energy behaviour, 

(77) 
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independent of its vector /axial- vector nature, its couplings, and the particles 
being scattered. This is the dominant contribution, even when y B is not small. 

Turning now to QCDC, we recall that the total cross-section is the weighted 
sum of parton-model pieces where there is no pt, and where there is a pt of 

2 _ \r\ 2 ^A. — n 2 ^ ~ x p^(^~ ~ z v) x v z v (jo\ 

whose relative importance is reduced by Q 2 /(Q 2 + 4p^). 

Studying the first-order cross-section (p0[), it is clear that only \M.2(fi, T2) | 2 
generates any 0-dependence, and that this is identical to that in the parton 
model with intrinsic p t given by (|78|) . One can therefore say that QCD does 
not play a direct part in determining the azimuthal correlations. Its role can 
be summarised as providing a large 'intrinsic' p t , and then diluting the resulting 
correlation. Furthermore, as mentioned above, this is dominated by the fact that 
the exchanged particle has spin-1. This also applies to the BGF part, although 
the dilution is not so great, because both contributions have an 'intrinsic' pt- 

We can calculate the size of the correlation, 

/ aa fda 3 cos (/) 

(cos0) = -jj^-, (79) 

(cos 20) = ^f^. (80) 

where the integrals are over whatever phase-space region is selected for the analy- 
sis. We examine two regions in particular, firstly the point where the correlation 
is maximised. Studying (^), we find that both (cos0) and (cos 20) are max- 
imised at x p = z p = 1/2. This corresponds to scattering in which both outgoing 
partons are at 90° to the boson-hadron axis in the Breit frame, with energy Q/2, 
so = 2p t = Q. Going to higher p t increases the size of the correlation gener- 
ated by \M.2(ri,r 2 )\ 2 , but also increases the dilution, so that the correlation gets 
smaller overall. At that point we have 

/ i \ -2y/l - y B (2 - y B ) 

K Wqodo.— - 9(l + (l-ite)a)+4(l- WB )+-Uife(2- Ifa )' 1 } 



9(1 + (1 - y B ) 2 ) + 4(1 - y B ) + AAy B {2 - y B ) ' 



(cos20) QCDCj _ = nM ■ M aL_ ^_ . __ _ (82 ) 



The BGF case is not quite so simple, because although the contributions to (cos 0) 
from the R\ and R2 terms in (|69|) are separately maximised at x p = z p = 1/2, 
they are equal and opposite there, so exactly cancel. The sum is maximised at 



x p = z p = (2 + ^8^3 — 12)/4 « 0.841, which corresponds to a configuration 
where one of the partons is at 90° to the boson-hadron axis, and the other is 
about 20° from the current direction, with p t « 0.159Q and y^s « 0.435Q- The 
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value at that point does not have a simple form, but is numerically similar to 
(cos 0) QCDC max - It would not be as easy to measure though, since different parts of 
phase-space give correlations with different signs, so integrating over some region 
of phase-space reduces (cos0) faster for BGF than for QCDC. For (cos 2(f)), the 
maximum is again at x p = z p = 1/2, with value 

{ms = i + (i-Ip+ B 4(i- fc r (83) 

independent of A, almost four times larger than (cos 20) QCDC max , and comparable 

to (cOS0) nr , nr 

i T I QCDC, max 

As already mentioned, as the transverse momentum is increased above Q/2, 
the correlation becomes weaker again. This is because the correlation generated 
by | A^tji; r 2) | 2 is governed by the direction of p 2 in the Breit frame, rather 
than its p t . Thus the correlation cannot get any stronger than it is when p 2 is 
at 90° to the current direction. However, the dilution by the azimuthally flat 
term increases with increasing s, so the correlation decreases overall. Since most 
experimental analyses use a fixed cutoff in p t irrespective of Q, it is primarily this 
dependence that determines the strong overall dependence on Q 2 . 

The other phase-space we consider is to not make any cuts at all. This is 
possible because both (cos0) and (cos 20) are infrared safe quantities, so if 
is measured semi-inclusively for all DIS events, the expectation values can be 
expanded as simple power series in a s [J We can then write the expectations 
symbolically as 

/ a\ Jda 3 cos0 

(cos0) = — — , (84) 

N ' 1 + / da 3 - Ca s v ' 

with an equivalent expression for (cos 20). Because the integrals are now over 
the whole phase-space, the integral in the denominator is divergent, but C, the 
first-order virtual correction to the lowest-order result, is also divergent such 
that after renormalisation and factorisation their sum is the finite first-order 
QCD correction to the total DIS cross-section. The integral in the numerator is 
convergent. Since it is first-order in a s , we can ignore the a s dependence of the 
denominator to first order, and we obtain simply 

(cos 0) = / da 3 cos 0. (85) 



+This is only strictly true if we define cj> m such a way that the expectation values are 
identical in the pt — > limit and at pt — 0, ie. if they are zero in the latter case. This can be 
achieved by adding the rule that if no particles in the event have any momentum transverse to 
the boson-hadron axis, a (f> value is chosen randomly and uniformly. Since experimentally such 
events never exist, this makes no practical difference to the result. 
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For QCDC we obtain 

/ * C v a s [i [i 1 9tK9-<S 

(cos (b) n „„„ = : — / ax„ / a^n 



QCDC 4vr 7, B p 7o p (l-x p )(l-^)x B /(x B ,g 2 ) 

., ., ., . .Av^ 2 - 1 cos # 2 sin 9 2 + 2l\/P - 1 sin 9 2 /0 ^ N 
x p (x 2 + x ± ) — — — , (86) 



which can be integrated over z p to give 




( cos 0>qcdc = — ' dx 



1 + Al + l 2 

which must be integrated numerically. We similarly obtain 

C F a s l 2 -l /-i x ^m,Q 



2\ 



(cos20) QCDC = — — — — / dx p x. 



,QCDC 4vr l + M + l 2 J XB v p x B f(x B ,Q 2 ) 
and 



cos</>)„^ = — ^— / dx n \/xJl — x 



^VF^T + 2/yF^l(2a: p - 1) 
1 + „4Z + Z 2 



(87) 



*9) 



Note that these are the extreme values of the expectations, assuming that the 
scattered parton direction could be perfectly identified. They would be reduced 
by a realistic method, such as using the hadron with largest Feynman-x F , or using 
all particles weighted by x F . The values of these expectations are shown in Fig. ^ 
in comparison with those from the intrinsic pt in the parton model. 

Although we have seen that the physical origin of the correlation is the same 
in QCD as in the parton model, the dependence on x B and Q 2 is quite difference, 
so the two contributions can be easily separated. 



5 The Drell-Yan Process 



The treatment of the Drell-Yan process again proceeds along similar lines. We 
start with the annihilation+gluon process, qq — > gV*, for which the matrix ele- 
ment is 



M 3 (X,p) = ig s vx{p 2 



(#2 ~ jQ 

—2p 2 -k 



-2p v k 



«a(PiW?)» ( 91 ) 
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10~ 2 10 _1 50 100 500 

x B Q a /GeV 2 

Figure 2: The average values of cos and cos 2<p from the parton model (dashed) 
and QCD (solid), for ep collisions at s = 10 5 GeV 2 with Q 2 = 10 2 GeV 2 (left) and 
x B = 0.01 (right). We use the MRS D-' distribution functions, A QCD = 230 MeV, 
with pure photon exchange, and (p t ) = y/ (p 2 ) = 500 MeV. The QCD curves are 
also broken down into the separate contributions from QCDC (dot-dashed) and 
BGF (dotted). 



where the momenta of the quark and antiquark are pi and p 2 . This can again be 
obtained from the QCD Compton case, by crossing 

P2 — > -P2- (92) 

The special vectors are 

(93) 
(94) 
(95) 
(96) 

We then obtain 



ri = 








r 2 = 


q-n 


= p 2 -k + p 1 


-Pi/xi 


T2 = 


P2/X2, 






fi = 


q-r 2 


= Pl -k + p 2 


-P2/X2 


\C + 


2 


8na s 





\C-\ 2 = 7 *™' ^A - (98) 



The cross-section in lowest order is 



da 2 = -dQ 2 dy / 9/1 (r ?1 )/ 5/2 (r ?2 )^|A^ 2 | 2 dr 2 , (99) 
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where y is the rapidity of the gauge boson, 771,2 = QI\Jse ±v are the energy- 
fractions of the quarks in the hadrons, f q /h are the corresponding distribution 
functions, and dT 2 is an element of the phase-space for the gauge boson decay. 
The corresponding first-order cross-section is 

11 C F a s dx 1 dx 2 d(p 

da 3 = -dsdy / g /i(fi)/$/ 2 (f2 



s » ^/ivw^/As^ 2 g 2 2n ( Xl + x 2 - 1) 3 (X! - l)(x 2 - 1) 2tt 
{x 2 1 \M 2 (r 1 ,f 2 )\ 2 + x 2 2 \M 2 (f 1 ,r 2 )\ 2 }dT 2 . (100) 

Just as in the e + e~ case, this can be evaluated exactly, without any knowledge 
of |-M 2 | 2 - The event can be boosted from the lab frame to the rest frame of 
the boson in such a way that the azimuth of the gluon around either pi or p 2 is 
identical to its lab-frame azimuth around the beam direction. However, these two 
different possibilities correspond to different configurations of the boson decay in 
the lab frame. \M. 2 {r\, f 2 )\ 2 clearly has no dependence on the azimuth around pi, 
and \A4 2 (ri, r 2 )\ 2 has no dependence on the azimuth around p 2 . This means that 
the boson decays can be correctly generated by generating a uniform azimuth 
around one or other of the quarks in the boson's rest frame, just as in e + e~ 
annihilation. 

It then only remains to relate the two cross-sections. Switching to using Q 2 
and t — — 2p 1 -k = —Q 2 (x 2 — 1) to parametrise the hadron-plane momenta, we 
obtain 

da 3 = \dsdydQ 2 di f q/1 ^)f m {i 2 ) ±. _1_ d ± 

{(Q 2 - u) 2 \M 2 ( ri ,f 2 )\ 2 + (Q 2 - i) 2 \M 2 (n, r 2 )\ 2 } dT 2 (101) 

and hence 

, / g /i(£i)/g/2(6) C F a s dsdi f 2 2 d<p 2 f 2 d(j)\ 
da z = t( Yfl \ o W - u ) da 2 — + {Q -t) da 2 — \ . 102 

fq/l{Vl)fq/2{V2) 2tT S 2 tu{ 2?T 2tT J 

The Monte Carlo algorithm is then 

1. Generate a parton model event according to the exact 0(1) matrix element. 

2. Generate s and i values according to Q 4 dsdi/s 2 iu = Q 4 dsdi/s 2 i(Q 2 — s — i). 

3. Construct momenta in the boson rest-frame corresponding to these values. 

4. Choose parton 1 with probability (Q 2 - uf / ((Q 2 - u) 2 + (Q 2 - t) 2 ) , and 
otherwise parton 2, and uniformly rotate about the chosen parton. 

5. Boost the momenta back to the lab frame such that the boson's rapidity is 
the same as it was in the parton-model event. 
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6. Calculate a weight factor 

w = ^ MM) { (1 _ M/g2)2 + (1 _ W) 2} (103) 

2?T J q /l{r}l)Jq/2{r}2) ^ } 

7. Keep the event with probability proportional to w. 

This correctly gives the properties of annihilation+gluon events, for example 
correlations between the gluon and lepton directions, with no knowledge of the 
matrix element that determines them, provided a generator of lowest-order events 
is available. 

Turning to the Compton process, qg — > qV*, we find that exactly the same 
treatment holds, and we obtain 

das = fl — \7t — \ o ^ W - t) d(T 2— + (Q -s) da 2 — \ , 104 

where t is this time —2p 2 -k. Equivalent formulae hold when the quark is replaced 
by an antiquark, or comes from the other hadron. 
The Monte Carlo algorithm is then 

1. Generate a parton model event according to the exact 0(1) matrix element. 

2. Generate s and t values according to —Q A dsdi/sH. 

3. Construct momenta in the boson rest-frame corresponding to these values. 

4. Choose parton 1 with probability (Q 2 - t) 2 / ((Q 2 - t) 2 + (Q 2 - s) 2 ) , and 
otherwise parton 2, and uniformly rotate about the chosen parton. 

5. Boost the momenta back to the lab frame such that the boson's rapidity is 
the same as it was in the parton-model event. 



6. Calculate a weight factor 

\a s / g /i(fi)/»/2(6) 



w 



t/Q 2 ) 2 + (l-s/Q 2 ) 2 }. (105) 



2vr f q /i(r]i)fg/2(V2) 

7. Keep the event with probability proportional to w. 

This again gives the correct event properties, with no knowledge of the matrix 
element that determines them, provided a generator of lowest-order events is 
available. 
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6 Summary 



By working in the simple gauge introduced by the CALKUL collaboration, it 
is possible to exactly factorise the first order corrections to the electroweak pro- 
duction, scattering and annihilation of quarks. This allows a transparent under- 
standing of event properties and correlations, and enables simple Monte Carlo 
prescriptions to be constructed. This was done in [jj for production processes, 
and in the present paper for scattering and annihilation processes. 

For production and annihilation the factorisation is complete, in the sense that 
the Monte Carlo algorithm does not need any details of the lowest-order cross- 
section to generate the correction. This is because it is able to use the angular 
configurations generated at lowest order, provided that all such configurations are 
generated. In the scattering case, the equivalent requirement is that all transverse 
boosts are generated, which is not a natural situation. Instead, the lowest order 
matrix element must be reevaluated for a new set of momenta. 

As an example of the additional insight that our approach can give, we dis- 
cussed azimuthal correlations between the leptonic and hadronic planes in DIS. It 
was seen that the correlation that arises in first-order QCD has the same physical 
origin as that in the parton model — simply the variation of the parton-lepton in- 
variant mass with azimuth, when there is a non-zero transverse momentum. The 
role of QCD is simply to provide that transverse momentum. By separating the 
cross-section into two components, our method is able to describe this transverse 
momentum in a gauge-invariant way. 

We finally note that since our results for first-order cross-sections are explicitly 
formulated as corrections to lowest-order cross-sections, they are ideally suited 
for making first-order matrix-element corrections to parton shower algorithms ||. 

References 

1. R. Kleiss, Phys. Lett. 180B (1986) 400. 

2. H. Georgi, H.D. Politzer, Phys. Rev. Lett. 40 (1978) 3. 

3. J. Chay, S.D. Ellis, W.J. Stirling, Phys. Rev. D45 (1992) 46. 

4. The CALKUL Collaboration, Nucl. Phys. B206 (1982) 53. 

5. R.N. Cahn, Phys. Lett. 78B (1978) 269. 
R.N. Cahn, Phys. Rev. D40 (1989) 3107. 

6. M.H. Seymour, LU TP 94-12, contributed to the 27th International Con- 
ference on High Energy Physics, Glasgow, U.K., 20-27 July 1994. 

M.H. Seymour, Matrix- Element Corrections to Angular- Ordered Parton 
Showers, in preparation. 



19 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/94 1 0244v 1 



